Astronomy & Astrophysics manuscript no. 14569rev © ESO 2010 

May 11,2010 



o 

o 

(N 



o 
u 

o 



(N 
> 

O 
(N 

in 

o 
o 



- 1—1 

X 



Letter to the Editor 

//^r5c/z^/ observations of FIR emission lines in briglitest cluster 

galaxies ^ 

A. C. EdgiS, J. B. R. OonlP, R. Mittai^l, S. W. AUeil^l, S. A. BaunP, H. BohringeiE], J. N. BregmaiP, M. N. BremeiEl, 

R CombeP, C. S. CrawforcP, M. DonahueffSl, E. EgamJUl, a. C. FabiaiJ^l, G. J. FerlancfEl, S. L. HameiO, N. A. 
HatcllE], W. JaffJ21, R. M. JohnstonJS, B. R. McNamariHl, C. R O'DeiEl, R PopessJ^, A. C. QuilleillS], r SalomiSl, 

C. L. SaraziiffS, g. M. VoiPSl, r. j. WilmaiffSl, and M. W. WisilSl 

' Institute for Computational Cosmology, Department of Physics, Durham University, Durham, DHl 3LE, UK 

^ Leiden Observatory, Leiden University, P.B. 9513, Leiden 2300 RA, The Netherlands 

^ Chester F. Carlson Center for Imaging Science, Rochester Institute of Technology, Rochester, NY 14623, USA 

'' Kavli Institute for Pai'ticle Astrophysics and Cosmology, Stanford University, 452 Lomita Mall, Stanford, CA 94305-4085, USA 

' Max-Planck-Institut fiir extraterrestrische Physik, 85748 Garching, Germany 

^ University of Michigan, Dept. of Astronomy, Ann Arbor, MI 48109, USA 

^ H H Wills Physics Laboratory Tyndall Avenue, Bristol BS8 ITL, UK 

^ Observatoire de Paris, LERMA, CNRS, 61 Av. de I'Observatoire, 75014 Paris, France 

' Institute of Astronomy, Madingley Rd., Cambridge, CB3 OHA, UK 

'" Michigan State University, Physics and Astronomy Dept., East Lansing, MI 48824-2320, USA 
" Steward Observatory, University of Arizona, 933 N. Cherry Avenue, Tucson, AZ 85721, USA 
'^ Department of Physics, University of Kentucky, Lexington KY 40506 USA 

'^ School of Physics and Astronomy, University of Nottingham, University Park, Nottingham NG7 2RD, UK 
''' Department of Physics & Astronomy, University of Waterloo, 200 University Avenue West, Waterloo, Ontario, Canada N2L 3G1 
'^ Department of Physics, Rochester Institute of Technology, 84 Lomb Memorial Drive, Rochester, NY 14623-5603, USA 
'^ Department of Physics and Astronomy, University of Rochester, Rochester, NY 14627, USA 
" Department of Astronomy University of Virginia, P.O. Box 400325, Charlottesville, VA 22904-4325, USA 
"* School of Physics, University of Melbourne, Victoria 3010, Australia 
" ASTRON, Netherlands Institute for Radio AstronomyP.O. Box 2, 7990 AA Dwingeloo, The Netherlands 

Received 30 March 2010/Accepted 

ABSTRACT 

The question of how much gas cools in the cores of clusters of galaxies has been the focus of many, multiwavelength studies in the 
past 30 years. In this letter we present the first detections of the strongest atomic cooling lines, [Cii], [Oi] and [Nii] in two strong 
cooling flow clusters, A1068 and A2597, using Herschel PACS. These spectra indicate that the substantial mass of cold molecular gas 
(> IO'Mq) known to be present in these systems is being irradiated by intense UV radiation, most probably from young stars. The 
line widths of these FIR lines indicate that they share dynamics similar but not identical to other ionised and molecular gas traced by 
optical, near-infrared and CO lines. The relative brightness of the FIR lines compared to CO and FIR luminosity is consistent with 
other star-forming galaxies indicating that the properties of the molecular gas clouds in cluster cores and the stars they form are not 
unusual. These results provide additional evidence for a reservoir of cold gas that is fed by the cooling of gas in the cores of the most 
compact clusters and provide important diagnostics of the temperature and density of the dense clouds this gas resides in. 

Key words. Galaxies: clusters: intracluster medium. Galaxies: clusters: elliptical and lenticular, cD 

1 . Introduction of gas cooling is related to energy injection into the intracluster 

medium by the action of jets and related AGN activity. 
The cooling process at the cores of galaxy clusters is highly 

complex: recent XMM-Newton and Chandra observations indi- ^he detection of substantial masses of molecular gas in the 
cate that the cooling rates are reduced by an order of magnitude ^^^^^ ^f ^^^ ^^^^ ^.^pi^jy ^^^jj^g ^j^^^^^^ ^^^^^^^ ^O lines 
below the simple cooling flow models at temperatures below (g^g^ 2OOI, Salome & Combes 2003) and warm H2 molecu- 
~ 2x lO^K (Peterson & Fabian 2006). These X-ray observations, j^ jj^^^ j^ ^^e NIR and MIR (Jaff^e & Bremer 1997, Egami et al. 
when linked with the detection of radio jet inflated bubbles in the 20O6) indicates that not all cooling is suppressed and this cooled 
cores of many of the strongest cooling flows (see McNamara & g^^ ^^^ ^^^^-^^^ ^^^ f^^j f^^ f^^^^^ ^q^ ^^^j^i^y ^^^^^ ^^^^^^^ 
Nulsen 2007 for a review), suggest that the strong suppression ^f molecular gas appear to correlate with the strength of optical 
lines from ionised gas (Crawford et al. 1999, Edge 2001) and the 



* Herschel is an ESA space observatory with science instruments dust continuum at MIR and sub-mm wavelengths (O'Dea et al. 
provided by European-led Principal Investigator consortia and with im- 2008). However, the excitation of these various emission lines 
portant participation from NASA. and the relative importance of energy input from star formation. 



A. C. Edge et al.: Herschel observations of FIR emission lines in brightest cluster galaxies 



AGN, cosmic rays and/or the intracluster medium is poorly con- 
strained (Ferland et al. 2009). 

One as yet unexplored diagnostic of the properties of the cold 
gas are the atomic cooling lines found in the FIR, [Cii], [Oi] and 
[Nil]. The unprecedented sensitivity of Herschel (Pilbratt et al. 
2010) to FIR line emission offers the opportunity to assess the 
ionisation and density of the colder gas for the first time with the 
[Cii] line and two principle [Oi] lines. The authors were awarded 
140 hours of time in an Open Time Key Project (PI Edge) to 
investigate the FIR line and continuum properties of a sample 
of 11 brightest cluster galaxies (BCGs) in well-studied cooling 
flow clusters selected on the basis of optical emission line and 
X-ray properties. The full goals of the project are to observe 
at least five atomic cooling lines for each object that cover a 
range in density and temperature behaviour and obtain a fully 
sampled FIR spectral energy distribution for systems where sig- 
nificant star formation is expected. In this paper we present the 
Photodetector Array Camera & Spectrometer (PACS, Poglitsch 
et al. 2010) spectroscopy for the two targets observed in the 
Science Demonstration Phase (SDP), Abell 1068 (z = 0.1386) 
and Abell 2597 (z = 0.0821). In a parallel paper (Edge et al. 
2010), we present the FIR photometry for these clusters. 

The two clusters observed have quite contrasting multiwave- 
length properties. A1068 is a strong MIR source (O'Dea et al 
2008) with a bright CO detection (Edge 2001) but a weak ra- 
dio source (McNamara et al 2004). A1068 lies just below the 
luminosity threshold of a ULIRG (10'^ L©) and exhibits some 
contribution from an AGN (Crawford et al. 1999, O'Dea et al. 
2008). On the other hand, A2597 is a relatively weak MIR source 
(Donahue et al. 2007) with a weak CO detection (Salome, priv. 
comm.) but a powerful radio source (Sarazin et al. 1995). The 
implied FIR luminosity of A2597 of 8.8x10'' Lq is a factor of 
around 30 below that of A1068 (3.5x10" Lq) and, in addi- 
tion, the fractional contribution from an AGN in the MIR is also 
lower. 

2. Observations 

We have observed the [Cii] and [Oi] lines at 157.74/im and 
63.18yLim for A 1068 and A2597 with the PACS spectrometer on 
Herschel. These are the primary cooling lines of the cold gas 
at a temperature T<40 K (Kaufman et al. 1999). In addition for 
A2597 we observed the [Nii], [Oi] and [Oiii] lines at \2\.9Q^iV(\, 
145.52/im and 88.36/im. These lines are used to constrain the 
excitation and temperature of this gas. Table 1 gives a summary 
of the observations. 

All spectral line observations were taken in PACS chopped 
line scan (standard faint line) mode with chopping-nodding. 
The simple pointed observations mode was used for all obser- 
vations. The data were reduced following the PACS data reduc- 
tion guide (PDRG) using the PACS Line Spectroscopy script 
for Point Source Chop/Nod Mode as presented by the PACS 
ICC team during the Herschel science demonstration phase data 
processing workshop at ESAC in december 2009. The reduc- 
tion was performed within the Herschel Interactive Processing 
Environment (HIPE) version 2.0.0 (Ott 2010), build RC 3. We 
have processed the data from level (raw channel data) to level 
2 (calibrated spectra) in a number of steps as outlined in the 
PDRG. Level to 1.0 processing removes the telescope specific 
structures from the data. The slopes of the raw channel data are 
fitted and removed. The signal is converted from data units to 
volts per second. Sky coordinate information is added and bad 
pixels and glitches are removed from data. The data is flatfielded 
and flux calibrated by applying the ground based nominal re- 



sponse function as recommended in the PACS spectroscopy per- 
formance and calibration (PSPC) document. This ground based 
response calibration is known to yield overestimated fluxes and 
following the PSPC we divide our fluxes by 1.3 and 1.1 in the 
blue and red bands. The accuracy of this flux calibration for the 
PACS spectrometer, at the time of writing, is about 50 percent 
within a given spectral band (PSPC). 

During the final stage of the reduction, level 1 .0 to 2.0, the 
data are spectrally and spatially rebinned into a 5x5xlambda 
cube. Using the standard 5x5 spatial rebinning each spatial pixel 
(spaxel) has a projected size of 9.4" x 9.4" on the sky. The spec- 
tral rebinning is performed using the recommended weak line 
density i.e. oversamp=l and upsamp=4. Values between 1 and 
10 were tried for the upsamp and oversamp parameters to test 
the robustness of the line profiles. We find that the line profiles 
do not change significantly for this range in values. 



3. Results 

The [Cii] 157/im and [Oi] 63jum lines are detected at a signal 
to noise greater than 30 for both A1068 and A2597. The much 
weaker [Nii] 122/im and [Oib] 145/im lines are detected at the 
3-5cr level for A2597. The [Oiii] 88;um line was not detected in 
A2597, an upper limit for this line is given in Table 2. 

The line spectra are fitted by a model consisting of; (i) a lin- 
ear function to determine the continuum flux, and (ii) a single 
gaussian function to determine the line flux. Continuum sub- 
tracted line spectra are shown for the central spaxel in Fig. 1 
for A1068 and Fig. 2 for A2597. The fitted line centers agree 
well with the redshift of CO in the BCG and the fitted FWHM 
line widths indicate gas with velocities of 300-500 km s ' . 

The [Cii] and [Nii] lines in the central spaxel of both ob- 
jects are well described by a single gaussian. However, the [Oi] 
63/im lines, where the PACS spectral resolution is best, have pro- 
files indicative of weak (2-3cr) deviations from a single gaussian 
function. The [Oi] line in A 1068 hints at a two-component struc- 
ture in the form of a narrow core component on top of a broad un- 
derlying component comparable to the CO(2-l) profile in Edge 
(2001). Both [Oi] lines observed in A2597 appear to have their 
dominant flux component at the systemic redshift of the BCG 
and a weaker component offset by about +250 km s"' which is 
also seen in the CO data (Salome, priv. comm). We attribute the 
shared structure of these atomic and moledular lines to gas kine- 
matics rather than self-absorption as the observed emission is 
from a large number of clouds that have much narrower intrinsic 
line width. 

The resolution of PACS at the observed wavelengths varies 
from about 5" for the [Oi] 63;um line to about 14" for the [Cii] 
157/im line. We have investigated line emission in all 25 spaxels 
of the PACS FoV. In all cases the line flux is dominated by the 
central spaxel. Summing up the flux in all 25 spectra and com- 
paring it to the flux in the central spaxel shows no evidence of 
excess line flux as compared to what is expected from a point 
source. In order to properly recover the full beam line fluxes we 
have applied point source corrections (appendix A of the PSPC 
document) to the central spaxel integrated line fluxes. The results 
are listed in Table 2. This spatial resolution matches the best sub- 
mm interferometry results for CO (Edge & Frayer 2003; Salome 
& Combes 2004) which implies that most of the emission is on 
scales < 5" so we believe our PACS line fluxes can be compared 
to literature values without large beam corrections. 
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Table 1. Log of Hers c he IP ACS observations 



Cluster 


Redshift 


Line 


Wavelength 
Oum) 


Obsid 


Bandwidth 
(kms"') 


Resolution 

(kms-') 


Beam Size 


A 1068 


0.1386 


[Cii] 


179.61 


1342186308 


1200 


201 


13.5"/33 kpc 






[Oi] 


71.94 


1342186307 


600 


55 


5.4"/13 kpc 


A2597 


0.0821 


[Cii] 


170.78 


1342187125 


1100 


218 


12.8720 kpc 






[Oi] 


68.41 


1342187124 


550 


68 


5.177.8 kpc 






[Nil] 


131.94 


1342188942 


1200 


281 


9.9"/15 kpc 






[Oib] 


157.56 


1342188704 


1200 


241 


11.8718 kpc 






[Om] 


95.61 


1342188703 




108 


7.2711 kpc 



Table 2. Spectral line results for HerschelPACS observations 



Cluster Redshift Line 



Integrated Line Flux 
(10-'*Wm-2) 



Velocity offset 
(kms-') 



measured FWHM 
(kms-') 



instrinsic FWHM 
(kms-') 



A 1068 


0.1386 


[Cii] 


104.7±1.8 


+25+55 


378+40 


320+55 






[Oi] 


64.8±0.2 


+25+50 


356+40 


352+55 


A2597 


0.0821 


[Cii] 


58.5+1.9 


+ 15+60 


463+40 


408+55 






[Oi] 


54.7+0.2 


+40+55 


411+40 


405+55 






[Nil] 


3.8+1.3 


+30+60 


578+90 


505+110 






[Oib] 


3.3+1.3 


-57+65 


484+90 


420+110 






[Om] 


<2.9 


- 


- 


- 



AIO68 CI 




179.5 IBO.O 

Wovelength (^m) 



A1068 01 




71.80 71.85 71.90 71.95 72.00 72.05 
Wovelength (fim) 



Fig. 1. Herschel PACS spectra of [Cii] and [Oi] in A1068 



4. Discussion 

The primary result from the SDP observations for this project is 
that the atomic cooling lines are present in both observed clus- 
ters. This first detection of these lines in cluster cores reinforces 
the importance of the cold gas in these environments. However, 
there are a number of questions that these detections raise. 

How do the properties of the FIR lines compare to lo- 
cal LIRGs/ULIRGs? There have been several studies of local 
galaxies with ISO and high redshift galaxies using ground- 
based instrument that cover [Cii] and [Oi] (Malhotra et al. 1997, 
Maiolino et al. 2005, Hailey-Dunsheath et al. 2010). These stud- 
ies show that the ratio of [Cii] to FIR luminosity is a function of 
luminosity with relatively less [Cii] emission for the most FIR lu- 
minous sources. Using the FIR data from Edge et al. (2010), we 
calculate the [Cii] to FIR luminosity ratios are lO"^"^ and lO-'"^ 
for A1068 and A2597, respectively. The [Cii]/FIR ratios of these 
two galaxies are comparable to those measured for galaxies of 
similar Lfir (see Fig. 2 of Maiolino et al. (2005)). In particu- 
lar, the [Cii] to FIR luminosity ratio is lower for the more FIR 
luminous of the two galaxies. The ratio of [Cii] to [Oi] shows 
less variation (1.62 and 1.07) and is again consistent with other 
comparable galaxies (Luhman et al. 2003). Also our CO(l-O) to 



FIR luminosity ratios of 10"^^ and 10"^^ for A1068 and A2597 
are consistent with star-forming local galaxies (Malhotra et al. 
1997). So, despite potential differences in excitation, pressure 
and metallicity, the relative intensity of the atomic and molecular 
lines to the FIR luminosity do not distinguish the BCGs studied 
here from other FIR bright galaxies. 



How do the dynamics of atomic and molecular lines com- 
pare ? The relative velocity width of the atomic lines compared 
to the CO and MIR H2 lines can provide important diagnos- 
tics for the dynamics and energetics of the various gas trac- 
ers. From the line width alone the resolution corrected line 
FWHM widths for the [Cii] and [Oi] lines are a;330 km s"' and 
ai400 km S-' for A 1068 and A2597 respectively. This compares 
to 243±13 km s-' (Edge 2001) and 292+45 km s-' (Salome, 
priv. comm.) for CO(2-l) for A 1068 and A2597. In each case 
the FIR lines are a factor of ~ 1 .35 broader This is not due to any 
instrumental broadening in the PACS instrument as the two lines 
sampled have similar intrinsic line width despite being observed 
at very different resolution. Instead, this difference is more likely 
to be related to the lines being emitted from different regions 
within the BCG or in shocks. However, this clearly needs to 
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A2597 Cll 




A2597 01 



170.0 170.2 170.4 170.6 170.8 171.0 171.2 171.4 
Wovelength (/^m) 

A2597 Nil 





68.3 68.4 

Wovelength (/xm) 

A2597 Olb 




131.5 132.0 132.5 

Wavelength (jtim) 



156.8 157.0 157.2 157.4 157.6 157.8 158.0 158.2 
Wavelength (^m) 



Fig. 2. Herschel PACS spectra of [Cii], [Oi] (63/zm), [Nn] and [Oib] (145yum) in A2597. 



be tested in more systems and through direct comparison of the 
[Cii] and [Oi] extent with that of CO. 

How do the FIR line ratios constrain the gas properties? The 
relative strength of the FIR lines can constrain several key prop- 
erties of the gas phase that dominates the emission. The main 
constraint we can determine directly from our current data is 
from the [Cii] 158jum and [Oi] 63 and 145yum Hnes for A2597. 
Kaufman et al. (1999) present photodissociation region (PDR) 
model predictions for the [Oi] 145/vm/63/vm and [Cii] 158//m 
to [Oi] 63/im line ratios. Combining these two constraints for 
our observed [Oi] 145jum/63jum ratio of 0.06+0.02 and [Oi] 
63//m to [Cii] 158//m ratio of 0.94+0.05, we estimate a density 
of 10^-^±"-5cm-3 and an incident FUV flux of Go of 150-1000 
Habing units. These values of Go imply intrinsic FUV luminosi- 
ties of a; 2-5x10'*-' erg s"' if the clouds subtend 3-5 kpc. This is 
comparable to the observed FUV luminosities of these galaxies 
once dust absorption is taken into account (O'Dea et al. 2004). 



5. Conclusions 

These initial results from Herschel indicate that atomic cooling 
lines are present in the brightest cluster galaxies in cooling flow 
clusters. The intensity and velocity width of these lines is con- 
sistent with all the other observed tracers of cold gas in these 
systems implying they originate from the same population of 
clouds. The only apparent exception to this in our current obser- 
vations is that the FIR lines appear to be systematically broader 
than the CO lines impling that the relative intensity of these 
lines varies with position within the BCG. The results that will 
come from our Open Time Key Project for 1 1 BCGs will expand 
greatly on those presented here with more lines and a greater 



range of BCG properties. Beyond this, the potential for Herschel 
to illuminate the properties of the cold gas that may fuel cold nu- 
clear accretion in more distant clusters and local groups is vast. 
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